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Abstract.  A  moving  tropical  cyclone  is  an  intense  localized  source  of  surface 
wind  stress  and  wind  stress  curl  that  produces  a  significant  response  in  the  ocean 
environment,  especially  in  the  ocean  thermal  structure,  the  upper  ocean  currents, 
and  the  sea  surface  elevation.  Such  a  response  has  been  well  identified  in  the 
open-ocean  region,  but  not  in  the  coastal  ocean  region.  In  this  study  we  use 
the  Princeton  Ocean  Model  with  20  km  horizontal  resolution  and  23  sigma  levels 
conforming  to  a  realistic  bottom  topography  to  identify  the  response  of  the  South 
China  Sea  to  Tropical  Cyclone  Ernie  1996.  Results  show  strong  Rimi1a.rit.iPR  in  the 
responses  between  open  ocean  and  coastal  regions,  including  near-surface  strong 
asymmetric  response  such  as  divergent  currents  with  near-inertial  oscillations, 
significant  sea  surface  temperature  cooling,  biase  to  the  right  of  the  storm  track,  sea 
surface  depressions  in  the  wake  of  the  storm,  and  subsurface  intense  upwelling  and 
cooling  at  the  base  of  the  mixed  layer  to  the  right  of  the  storm  track.  The  unique 
features  of  the  SCS  response  to  Ernie  are  also  discussed. 


1.  Introduction 

Observational  and  modeling  studies  on  oceanic  re¬ 
sponse  (thermal,  circulation,  and  surface  elevation)  to 
tropical  cyclones  have  been  focused  on  open  ocean  re¬ 
gions  (ORs)  but  not  on  coastal  regions  (CRs),  espe¬ 
cially  the  South  China  Sea  (SCS).  The  ocean  thermal 
response  identified  in  ORs  is  the  sea  surface  cooling  that 
is  very  important  to  the  interaction  between  tropical  cy¬ 
clones  and  the  ocean  [Emanuel,  1988;  Chang,  1985;  Chu, 
1993;  Chu  et  al,  1990;  Chu  and  Garwood,  1991].  Obser¬ 
vational  evidence  clearly  show  that  sea  surface  temper¬ 
ature  (SST)  usually  decreases  by  several  degrees  Celsius 
because  of  the  storm  forcing.  For  example,  Fisher  [1958] 
documented  a  3°C  cooling  by  analyzing  ship  reports 
after  the  passage  of  Tropical  Cyclones  Connie  and  Di¬ 
ana  in  1955.  Jordan  [1964]  showed  that  maximum  SST 
decreases  occur  on  the  right  side  of  the  storm  track. 
A  warm  layer  occurred  near  the  surface  at  the  radius 
of  maximum  winds  (J?m),  then  deepened  with  increas¬ 
ing  distance  away  from  the  storm  center.  Black  [1983] 
also  observed  this  pattern  in  data  from  Typhoon  Ella  in 
1978.  He  reports  a  negative  annulus  of  wind  stress  curl 
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may  produce  the  warming  at  the  outer  regions  of  the 
storm.  Pudov  et  al.  [1978]  suggested  that  downwelling 
is  the  mechanism  that  produces  this  warm  layer. 

The  ocean  circulation  response  identified  in  ORs  is 
a  significant  right-hand  bias  in  the  upper  ocean  cur¬ 
rent  amplitude  with  more  than  1  m  s'1  in  the  mixed 
layer.  Price  [1981]  attributes  this  right-hand  bias  in  the 
current  amplitude  to  an  inherent  asymmetry  in  the  cou¬ 
pling  between  the  wind  stress  of  the  moving  cyclone  and 
the  wind-driven  mixed  layer  currents.  This  mixed  layer 
current  structure  produces  a  unique  convergence  and 
divergence  pattern  in  the  wake  of  the  storm.  Thermo- 
cline  currents  are  opposite  to  those  in  the  mixed  layer. 
Surface  divergence  and  upwelling  produces  a  subsur¬ 
face  pressure  anomaly  [Price  et  al.,  1994],  Subsurface 
current  converges  inward  toward  this  pressure  anomaly, 
producing  very  strong  current  shear  between  the  mixed 
layer  and  the  thermocline.  The  inertially  rotating  wind- 
driven  currents  are  highly  divergent,  which  forces  an  os¬ 
cillating  vertical  velocity  at  the  base  of  the  mixed  layer 
with  a  near-inertial  period.  Ocean  sensors  recorded 
the  first  velocity  structure  measurements  in  the  Gulf  of 
Mexico  during  the  passage  of  Hurricane  Eloise  [  Withee 
and  Johnson,  1988.]  This  data  set  clearly  shows  the 
generation  and  rapid  decay  of  these  near-inertial  mo¬ 
tions  at  53  m,  roughly  the  mixed  layer  depth,  before 
and  after  the  passage  of  the  storm.  Shay  and  Elsberry 


13,991 


Report  Documentation  Page 


Form  Approved 
OMB  No.  0704-0188 


Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 


1.  REPORT  DATE 

15  JUN  2000 

2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2000  to  00-00-2000 

4.  TITLE  AND  SUBTITLE 

Response  of  the  South  China  Sea  to  Tropical  Cyclone  Ernie  1996 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Postgraduate  School, Monterey  ,CA, 93943 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

1 

1 

16.  SECURITY  CLASSIFICATION  OF: 


a.  REPORT 

unclassified 


b.  ABSTRACT 

unclassified 


c.  THIS  PAGE 

unclassified 


17.  LIMITATION  OF 
ABSTRACT 

Same  as 
Report  (SAR) 


18.  NUMBER 
OF  PAGES 

19 


19a.  NAME  OF 
RESPONSIBLE  PERSON 


Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


°m  ET  AL:  RESP0NSE  OF  THE  S°™  «*  ™  ™or.cAL  cvctoNE  ERNIE  1996 

rino>7l  ... 


?,8?!  e,Xjmmed  a  more  extensive  set  of  observations 
collected  during  the  passage  of  Hurricane  Frederic  1979 
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topography  of  the  South  China  Sea  (SCS) 
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Tropical  Storm  Ernie  (Nov  1996)  Best  Track 
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Figure  2.  Track  of  Tropical  Cyclone  Ernie  1996. 


the  system  to  begin  moving  to  the  west-southwest  back 
into  the  SCS.  On  November  13,  Ernie  moved  back  to 
the  original  position  of  November  9. 

2.3.  West-Southwestward  Movement  Stage 
(November  13-18) 

Ernie  continued  moving  to  the  west-southwest  over 
the  next  4  days,  passing  over  the  southern  tip  of  Viet¬ 
nam,  and  slowly  weakening  to  15  m  s-1.  Ernie  tracked 
into  the  Gulf  of  Thailand  and  then  finally  into  the  east¬ 
ern  Bay  of  Bengal,  where  it  dissipated  over  water  on 
November  18.  The  storm  moved  southwesterly  at  5  m 
s_1.  The  storm  maintained  its  size  during  the  looping 
stage  (Rm  —  72  km);  however,  the  storm  strength  was 
reduced  to  moderate,  given  the  JTWC-analyzed  wind 
speeds  of  15-20  ms-1. 


where  r  is  the  horizontal  distance  to  the  storm  cen¬ 
ter,  (uc,ve)  are  the  radial  and  tangential  velocity  com¬ 
ponents,  7  is  the  inflow  angle  of  the  wind  as  it  spi¬ 
rals  into  the  center  of  the  cyclone,  a  is  a  scaling  factor 
(a  =  r( Rm)  that  makes  ue(r)  continuous  at  the  tropical 
cyclone  center,  and  X  is  a  positive  constant  <  1.  Con¬ 
sidering  the  beta  effect  propagation,  Carr  and  Elsberry 
[1997]  proposed  using  X  ~  0.4. 

To  prevent  a  large  discontinuity  at  the  periphery  of 
the  storm,  however,  we  blend  the  tropical  wind  field 
into  the  interpolated  5  day  mean  NASA  scatterometer 
(NSCAT)  wind  fields  on  November  5-9,  10-14,  and  15- 
19.  Adjustments  were  made  near  the  storm’s  periphery 
to  smooth  the  transition  between  the  background  and 
storm  wind  fields.  This  was  done  by  a  weighted  average 
inside  and  outside  of  the  tropical  cyclone, 


3.  A  Tropical  Cyclone  Wind  Profile 
Model 

3.1.  Model  Description 

TCWPM,  a  recently  developed  tropical  cyclone  wind 
profile  model  (TCWPM)  by  Carr  and  Elsberry  [1997], 
is  used  to  establish  a  high-resolution  surface  wind  field 
for  Tropical  Cyclone  Ernie  1996.  Here  let  Ro  and  Rm 
represent  the  radii  of  zero  and  maximum  tangential  ve¬ 
locities  inside  a  tropical  cyclone.  On  the  base  of  the 
angular  momentum  balance,  Carr  and  Elsberry  [1997] 
proposed  a  model  to  compute  the  wind  vector  relative 
to  the  center  of  the  tropical  cyclone, 


(1) 


V  =  (1  —  e)(Vc  +  Vt)  +  eVjj,  (2) 


where  V  kg  is  the  background  wind  field,  Vj  is  the  storm 
translation  velocity,  and  the  weight  e  is  computed  by 


e  — 


l+c4 


r 

0.9Ro 


(3) 


The  effects  of  these  adjustments  are  to  increase  the 
weighting  on  the  storm  wind  gradually  and  decrease 
the  weighting  on  the  background  winds  as  the  radial 
distance  to  the  storm  center  decreases.  Both  effects 
provide  for  a  smoother  transition  between  the  storm 
wind  field  and  the  background  wind  field.  The  daily 
evolution  of  the  blended  surface  wind  fields  is  shown  in 
Figure  3. 


3.2.  Wind  Model  Verification 

We  use  the  observational  NSCAT  data  to  verify  the 
model  wind  fields.  While  the  comparisons  were  per- 


uc(r)  =  tan(7)vc(r), 
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Figure  3.  Tropical  Cyclone  Wind  Profile  Model  (TCWPM)  generated  surface  wind  field  over 
SCS  on  November  6-17,  1996,  representing  Ernie. 


formed  for  the  history  of  Tropical  Storm  Ernie,  modeled 
wind  vector  field  (Figure  3)  comparisons  with  observa¬ 
tional  NSCAT  winds  (Figure  4)  were  only  conducted 
for  the  SCS  domain  98°-121°E  and  3®S-25°N  at  the 
observational  points.  Comparing  Figure  3  to  Figure 
4,  Ernie  was  well  represented  in  both  TCWPM  and 
NSCAT  winds.  However,  the  difference  between  the  two 
wind  fields  is  not  negligible.  Figure  5  shows  the  tempo¬ 
ral  variation  of  the  root-mean-square  (rms)  wind  vector 
difference.  The  largest  rms  errors  were  noted  during 
the  November  11-13  period  when  the  storm  was  mov¬ 
ing  over  Luzon,  where  significant  modifications  to  the 
wind  field  occurred.  The  wind  model  used  in  this  study 
does  not  account  for  land  interaction,  however,  and  this 


aspect  of  the  model  appears  to  be  causing  these  large 
errors. 

This  wind  field  error  analysis  suggests  that  with  the 
exception  of  the  period  that  Ernie  was  moving  over 
Luzon,  the  wind  model  was  successful  in  simulating 
the  observable  characteristics  of  the  wind  field.  How¬ 
ever,  there  is  uncertainty  in  this  analysis  as  the  NSCAT 
instrument  cannot  observe  the  entire  SCS  during  one 
pass. 

4.  The  Princeton  Ocean  Model 

4.1.  Model  Description 

The  POM  is  a  time-dependent,  primitive  equation 
circulation  model  on  a  three-dimensional  grid  that  in- 
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Figure  3.  (continued) 


eludes  realistic  topography  and  a  free  surface  [Blumberg 
and  Mellon ,  1987] .  The  horizontal  spacing  is  0.179° 
by  0.175°  (“20  km  resolution),  and  there  are  23  ver¬ 
tical  sigma  coordinate  levels.  The  model  domain  is 
from  3.06°S  to  25.07°N,  and  98.84°  to  121.16°E,  which 
encompasses  the  SCS  and  the  Gulf  of  Thailand,  and 
uses  realistic  bathymetry  data  from  the  Naval  Oceano¬ 
graphic  Office  Digital  Bathymetric  Data  Base  5  min  by 
5  min  resolution  (DBDB5).  Consequently,  the  model 
contains  T25  x  162x  23  horizontally  fixed  grid  points. 
The  horizontal  diffusivities  are  modeled  using  the  Sma- 
gorinsky  [1963]  form  with  the  coefficient  chosen  to  be  0.2 
for  this  application.  The  bottom  stress  rb  is  assumed 
to  follow  a  quadratic  law  [Blumberg  and  Mellor,  1987] 
in  our  model. 


4.2.  Atmospheric  Forcing 

tnT^"tmrheric  forcing  for  the  SCS  application  of 
the  POM  includes  mechanical  and  thermohaline  forcing. 
The  wind  forcing  is  depicted  by 

,,  ,du  dv. 

M  dz  '  dz^z=0  =  (r°x'  r°v)  (4) 

where  (u,  v)  and  (r0jr,r0y)  are  the  two  components  of 
the  water  velocity  and  wind  stress  vectors,  respectively. 
The  restoring  type  surface  thermohaline  forcing  is  used 
in  this  study  is 


—  C'C^OBS  -  9) 


(5) 
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Figure  4.  NASA  Scatterometer  (NSCAT)  surface  wind  field  over  SCS  on  November  6-17,  1996, 
representing  Ernie. 


=  C(5'obs  -  S’),  (6) 

where  $obs  and  -Sobs  are  the  observed  potential  tem¬ 
perature  and  salinity,  respectively.  The  relaxation  co¬ 
efficient  C  is  the  reciprocal  of  the  restoring  time  period 
for  a  unit  volume  of  water.  The  relaxation  coefficient 
C  is  taken  to  be  0.7  m  d-1,  which  is  equivalent  to  a 
relaxation  time  of  43  days  for  an  upper  layer  30  m  thick 
\Chu  et  al.,  1998b,  1999b,  c].  The  net  effect  is  to  pre¬ 
vent  any  deviation  from  climatology  and  ensure  that 
the  SCS  acts  as  a  heat  source. 


4.3.  Lateral  Boundary  Forcing 

Rigid  lateral  boundaries,  i.e.,  the  modeled  ocean  bor¬ 
dered  by  land,  were  defined  using  a  free  slip  condi¬ 
tion  for  velocity  and  a  zero-gradient  condition  for  tem¬ 
perature  and  salinity.  No  advective  or  diffusive  heat, 
salt,  or  velocity  fluxes  occur  through  these  boundaries. 
Open  boundaries  are  where  the  numerical  grid  ends 
but  the  fluid  motion  is  unrestricted.  Uncertainty  at 
open  boundaries  makes  marginal  sea  modeling  difficult. 
Three  approaches,  local-type,  inverse-type,  and  nested 
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Figure  4.  (continued) 


basin/coastal  modeling,  are  available  for  determining 
open  boundary  conditions.  Here  we  take  the  local- 
type  approach,  i.e.,  use  the  radiative  boundary  condi¬ 
tion  with  specified  volume  transport.  When  the  water 
flows  into  the  model  domain,  temperature  and  salinity 
at  the  open  boundary  are  likewise  prescribed  from  the 
climatological  data  [Levitus,  1982],  When  water  flows 
out  of  the  domain,  the  radiation  condition  was  applied, 

dt^9'  ^  +  Un  dn^9'  ^  =  0  (7) 

where  the  subscript  n  is  the  direction  normal  to  the 
boundary. 


Notice  that  flows  through  SCS  straits  are  quite  un¬ 
certain.  Choosing  one  among  various  estimations  is  not 
an  easy  job.  Despite  being  old,  Wyrtki’s  [1961]  data 
provide  a  balanced  estimation  of  volume  transports  for 
the  Luzon  Strait,  Taiwan  Strait,  and  Gasper/Karimata 
Strait  with  seasonal  variations  (Table  1).  Since  there 
are  no  reliable  estimations  at  the  Balabac  Channel,  the 
Mindoro  Strait,  and  the  Strait  of  Malacca,  we  assumed 
zero  transport  there.  Such  a  treatment,  especially  dur¬ 
ing  the  tropical  storm  forcing,  may  distort  the  solution. 
Applying  fixed  inflow  boundary  conditions  leads  to  a 
dynamical  inconsistency  and  the  development  of  spuri¬ 
ous  Poincare  and  Kelvin  waves  in  that  vicinity. 
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Figure  5.  Temporally  varying  rms  difference  (ms-1)  between  TCWPM  and  NSCAT  wind  vector 
components  (a)  u  and  (b)  v. 


4.4.  Experiment  Design 

The  model  year  consists  of  360  days  (30  days  per 
month);  day  361  corresponds  to  1  January.  The  inte¬ 
gration  was  divided  into  preexperimental  and  experi¬ 
mental  stages.  During  the  preexperimental  stage  the 
model  was  integrated  for  34  months  and  3  days  from 
an  initial  at  rest  state  with  three-dimensional  climato¬ 
logical  January  temperature  and  salinity  fields  [levitus, 
1982]  forced  by  the  climatological  monthly  mean  wind 
stress  [Hellerman  and  Rosenstein,  1983].  The  final  state 
from  the  preexperimental  run  was  taken  as  the  SCS 
condition  for  November  4,  1996,  defined  as  the  initial 
condition  for  the  experimental  stage.  During  the  exper¬ 
imental  stage,  the  model  was  forced  by  the  embedded 
TCWPM/NSCAT  wind  field  to  simulate  the  oceanic 
response  to  Tropical  Cyclone  Ernie,  for  18  days.  This 
wind  model  attempted  to  simulate  realistically  not  only 
the  movement  of  the  storm  but  also  its  size  and  intensity 
variability.  On  the  eighteenth  day  the  storm  passed  out 
of  the  model  domain  into  the  Bay  of  Bengal,  and  the 
model  run  was  stopped.  Three-dimensional  u,v,  and  w , 
temperature  and  salinity  fields,  two-dimensional,  depth- 
averaged  u  and  v  velocity  fields,  and  sea  surface  eleva¬ 
tion  fields  were  output  by  the  model. 

5.  Numerical  Simulation 

The  major  purpose  of  this  study  is  to  identify  whether 
the  following  strong  asymmetric  ocean  responses  to  a 
tropical  storm  still  exist  in  the  CR  regions:  (1)  surface 


and  subsurface  current  reversal,  (2)  near-inertial  asym¬ 
metric  currents,  (3)  divergence-convergence-induced  ver¬ 
tical  motion,  (4)  significant  sea  surface  depression  and 
cooling,  biased  to  the  right  of  the  storm  track,  and  (5) 
subsurface  intense  cooling  at  the  base  of  the  mixed  layer 
to  the  right  of  the  storm  track. 

5.1.  General  Features 

5.1.1.  Surface  and  subsurface  current  rever¬ 
sal.  Oceanic  current  responses  are  pronounced  dur¬ 
ing  Ernie’s  movement  over  the  SCS  on  November  7-18. 
Outflow  from  the  storm  produced  very  strong  diver¬ 
gent  upper  layer  currents  (Figure  6).  This  is  clearly 
seen  from  the  cross  section  of  the  v  velocity  field  taken 
along  12.32I>N  on  November  16  (Figure  7a).  Velocities 
were  20-200  cm  s-1  in  the  surface  layer.  The  model 
produced  very  intense  vertical  current  shear  across  the 
base  of  the  mixed  layer,  with  40  cm  s-1  northward 
flow  in  the  mixed  layer  and  40  cm  s'1  southward  flow 
in  the  thermocline  (Figure  7a).  Shay  et  al.  [1992] 
found  a  similar  pattern  in  current  profiler  measure¬ 
ments  from  Hurricane  Gilbert.  This  180°  phase  rever¬ 
sal  between  the  mixed  layer  and  thermocline  indicates 
that  wind  stress  generated  the  near-surface  currents, 
whereas  pressure  gradient  effects  generated  the  ther¬ 
mocline  currents  [Price,  1983]. 

To  assess  the  pressure  gradient  errors  over  steep  bot¬ 
tom  slope  [Haney,  1991],  we  present  the  cross  section 
of  the  v  velocity  field  taken  along  12.32°N  on  Novem¬ 
ber  4  (Figure  7b),  the  day  Ernie  was  generated  (Figure 


Table  1.  Bimonthly  Variation  of  Volume  Transport  at  the  Lateral  Open  Boundaries 


Feb. 

April 

June 

Aug. 

Oct. 

Dec. 

Gaspar-Karimata  Straits 

4.4 

0.0 

-4.0 

-3.0 

1.0 

4  3 

Luzon  Strait 

-3.5 

0.0 

3.0 

2.5 

-0.6 

-3.4 

Taiwan  Strait 

-0.9 

0.0 

1.0 

0.5 

-0.4 

-0.9 

The  positive/negative  values  mean  outflow/inflow  were  taken  from  Wyrtki  [1961],  Volume  transport  is  in  Sverdrups. 
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Figure  6.  Model-generated  surface  currents  on  November  8-17,  1996.  Here  open  circles  indicate 
the  location  of  Ernie. 


2).  The  pressure  gradient  errors  generated  by  the  steep 
topography  should  be  comparable  on  November  4  and 
16.  However,  we  only  see  the  western  boundary  cur¬ 
rents  on  November  4  but  not  the  flow  reversal  between 
the  mixed  layer  and  the  thermocline.  This  shows  that 
the  surface  and  subsurface  current  reversals  are  caused 
by  the  winds  from  Ernie. 

5.1.2.  Near-inertial  currents.  Storm-generated 
near-inertial  currents  were  previously  reported  for  the 
OR  regions  [e.g.,  Shay  et  al.,  1992],  Here  we  present 
the  time  series  of  current  components  (u  and  n)  at  cen¬ 
tral  SCS  (13°N,  119. 5°E)  for  two  depths:  30  and  150  m 
(Figure  8).  Velocities  at  30  m  depth  represent  the  up¬ 
per  ocean  currents,  and  velocities  at  150  m  denote  the 
thermocline  currents.  Velocities  at  both  depths  reveal 
nearly  2  days  of  oscillation  with  the  following  features: 


(1)  very  weak  oscillation  as  Ernie  approached  the  SCS 
(before  November  8),  (2)  strong  oscillation  after  Ernie 
entered  the  SCS,  (3)  flow  reversal  (180°  out-of-phase 
between  the  upper  ocean  currents  and  the  thermocline 
currents),  and  (4)  larger  amplitudes  of  u  and  v  time  se¬ 
ries  at  30  m  depth  (upper  ocean)  than  at  150  m  depth 
(thermocline).  The  power  density  spectrum  versus  the 
period  of  the  near-surface  u  and  v  velocities  (Figure  8) 
shows  a  peak  located  at  2. 1  days,  which  is  very  close 
to  the  local  inertial  period  at  13° N  of  2.2  days.  This 
implies  a  blue  shift  in  the  frequency  of  the  near-inertial 
motions  consistent  with  the  theoretical  argument. 

5.1.3.  Asymmetric  currents.  Another  striking 
surface  feature  produced  by  the  storm  during  its  move¬ 
ment  was  the  strong  “fan-shaped”  pattern  of  convergent 
and  divergent  currents  to  the  right  of  the  storm  track 
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Figure  6.  (continued) 


(Figure  6).  When  the  surface  currents  are  onshore  to 
the  east  and  southeast,  coastal  downwelling  occurred, 
and  a  strong  northwesterly  subsurface  jet  was  produced. 
When  the  currents  reversed  and  became  offshore  to  the 
west  and  northwest,  coastal  upwelling  occurred.  The 
subsurface  jet  reversed  to  the  southeast  and  weakened 
slightly.  The  alternating  convergence  and  divergence 
patterns  in  the  surface  current  fields  are  to  the  right 
of  the  storm  track.  The  450  km  half  wavelength  of 
this  current  pattern  agrees  quite  well  with  a  calculated 
theoretical  value  of  the  inertial  half  wavelength. 

5.1.4.  Surface  cooling  to  the  right  of  the  storm 
track.  Similar  to  hurricane  behavior  in  the  OR  re¬ 
gion,  Ernie  generated  surface  cooling  to  the  right  of  the 
storm  track.  SST  decreases  of  0.25°C  were  noted  to 
the  right  of  the  storm  track  (Figure  9),  which  is  con¬ 
sistent  with  previous  observations  of  tropical  cyclone 
SST  cooling  [Jordan,  1964;  Hazelworih,  1968;  Black, 
1983].  The  model  produced  a  crescent-shaped  pattern 
of  SST  decrease  to  the  right  the  storm  track.  This 
pattern  is  consistent  with  earlier  studies  by  Pudov  et 
al.  [1978]  and  Black  [1983].  Maximum  cooling  was 
1.5°C  at  80  km  from  the  track,  '1.5  During 

the  west-southwestward  movement  stage  (November  IS¬ 
IS),  Ernie  gradually  lost  intensity  as  its  wind  speed  de¬ 
creased  from  18  to  13  m  s-1.  Because  of  the  storm’s 
weak  intensity,  surface  responses  in  the  wake  of  the 
storm  were  weak.  Cooling  was  only  0.25°C  and  cov¬ 
ered  a  broad  area  to  the  right  of  the  storm  track. 


5.1.5.  Surface  depression  to  the  right  of  the 
storm  track.  Usually,  hurricane-generated  surface  de¬ 
pression  in  the  OR  region  is  confined  to  the  storm  track 
[Shay  et  al.,  1990].  Similarly,  Ernie-generated  surface 
depression  was  also  confined  to  the  storm  track.  Along 
the  coast  the  model  generated  an  oblong  sea  surface  de¬ 
pression  of  5-10  cm  by  November  8  (Figure  10).  Strong 
surface  divergence  and  coastal  effects  caused  by  storm 
wind  may  produce  this  depression.  By  00:00  LT  on 
November  9,  Ernie  was  located  near  16°N,  118°E.  The 
model  simulated  a  sea  surface  trough,  with  maximum 
depressions  of  10  cm,  confined  to  the  track,  in  the  wake 
of  the  storm  (Figure  10).  The  trough  extends  100-150 
km  on  either  side  of  the  storm  track,  2-3  R^.  During 
the  west-southwestward  movement  and  weakening  stage 
(November  13-18)  the  sea  surface  depression  was  only 
5  cm. 

5.1  6.  Subsurface  cooling.  The  model  fields  in¬ 
dicated  subsurface  cooling  in  the  mixed  layer  from  1 10° 
to  114° E  in  a  cross  section  at  12.3°N  on  November  16 
(Figure  11a).  Strong  surface  divergence  caused  maxi¬ 
mum  cooling  of  more  than  1.5°C  in  this  region  at  50-80 
m.  Along  the  coast,  warming  of  0.5°C  occurred  down 
to  depths  in  excess  of  200  m  because  of  intense  down- 
welling  produced  by  offshore  current  flow  near  13°  N. 
Twenty-four  hours  later,  these  patterns  were  reversed. 
Temperatures  within  a  narrow  1°  wide  wedge  of  the 
broad  cool  anomaly  warmed  significantly.  Temperature 
increases  of  up  to  2.5°C  occurred  (Figure  lib)  down  to 
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and  (b)  November  4  and  (c)  T  (degree  Celsius)  on  November  16,  1996. 


175  m.  Along  the  coast,  upwelling  caused  temperatures 
to  cool  by  PC.  An  analysis  of  the  surface  current  pat¬ 
tern  indicates  that  the  forward  propagation  of  the  near- 
inertial  surface  current  generated  by  the  storm  caused 
this  reversal  of  the  subsurface  temperature  anomalies. 
Near  the  coast  the  convergence  and  onshore  flow  re¬ 
placed  upper  ocean  divergence  and  offshore  flow. 

5.2.  Unique  Features 

The  model  also  simulated  several  unique  features,  ap¬ 
parently  caused  by  Kuroshio  intrusion,  coastal  interac¬ 
tions  under  storm  forcing,  and  quasi-stationary  storm 
forcing. 

5.2.1.  Convergence  between  coastal  currents 
and  intruded  Kuroshio.  Observational  studies  show 
an  evident  intrusion  of  Kuroshio  water  into  the  SCS 
through  Luzon  Strait  during  the  winter  monsoon  sea¬ 


son  [Wyrtki,  1961;  Shaw,  1991],  The  influence  of  Ernie 
on  the  Kuroshio  intrusion  was  simulated  by  the  POM 
as  a  SST  increase  (1°C)  across  the  Luzon  Strait  from 
southwest  Taiwan  to  northwest  Luzon.  This  region  ap¬ 
pears  to  be  a  convergence  zone  between  a  northward 
flowing  coastal  current  along  the  western  coast  of  Lu¬ 
zon  and  the  inflow  of  a  branch  of  the  warm  Kurishio 
Current  through  the  Luzon  Strait  (Figure  6).  Such  a 
convergence  produced  downwelling  and  raised  SST  be¬ 
tween  P  and  1.5*0  (Figure  9). 

5.2.2.  Coastal  jets  generated  by  storm  forcing. 
Coastal  interactions  under  storm  forcing  appear  to  have 
produced  a  subsurface  jet  ('40  cm  s-1)  flowing  north¬ 
ward  and  parallel  to  the  coast  of  Luzon  at  50-100  m  on 
November  9  (Figure  12).  A  time  series  of  the  u  and 
v  velocities  at  13°N,  119. 5°E  (Figure  8)  suggests  that 
this  feature  is  related  to  the  oscillating  surface  currents. 
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U  Time  Series  13.02N  119.5E 


Figure  8.  Time  series  of  (a)  u  and  (b)  v  and  (c) 
119.5°E  from  November  2  to  18,  1996. 


power  density  spectra  of  u  and  v  at  13°N, 


When  the  surface  currents  are  onshore  to  the  east  and 
southeast,  coastal  downwelling  occurred,  and  a  strong 
northwesterly  subsurface  jet  was  produced.  When  the 
currents  reversed  and  became  offshore  to  the  west  and 
northwest,  coastal  upwelling  occurred.  The  subsurface 
jet  reversed  to  the  southeast  and  weakened  slightly.  A 
divergent  surface  outflow  caused  sea  surface  rises  along 
the  coast  of  up  to  10  cm  extending  along  the  entire  west 
coast  of  Luzon  (Figure  10). 

5.2.3.  Concentric  response  to  looping  around 
the  storm.  During  Ernie’s  looping  around  in  the 
northern  SCS  on  November  9-13  (Figure  2),  SST  cool¬ 
ing  and  sea  surface  setdown  were  concentrated  along 
the  track  (Figures  9  and  10).  This  may  be  attributed  to 
interaction  with  the  near-inertial  oscillations  produced 
by  Ernie  during  its  first  pass  through  the  region.  As 
the  storm  initially  moved  through  the  region,  it  set  up 
near-inertial  oscillations  in  the  velocity  and  tempera¬ 
ture  fields.  Furthermore,  SST  and  sea  surface  elevation 
anomaly  patterns,  which  had  been  concentrated  in  the 
storm’s  wake,  quickly  became  nearly  concentric  with 
the  storm’s  position  as  it  went  quasi-stationary.  Maxi¬ 
mum  SST  cooling,  1°C,  and  sea  surface  depression,  10 
cm,  occurred  near  the  storm. 


6.  POM  Model  Verification 

6.1.  SST 

The  multichannel  sea  surface  temperature  (MCSST) 
data  obtained  from  the  polar  orbiting  Defense  Meteoro¬ 
logical  Satellite  Program  (DMSP)  is  used  to  verify  the 
POM-modeled  SST  fields.  Because  of  extensive  cloud 
cover  associated  with  Ernie,  MCSST  observations  are 
limited  to  cloud-free  regions  with  the  maximum  number 
of  data  points,  1180,  on  November  16  and  the  minimum 
number  of  data  points,  256,  on  November  9  (Figure 
13).  The  POM-simulated  SST  fields  agree  well  quali¬ 
tatively  with  the  DMSP  MCSST  fields  (Plate  1).  For 
a  quantitative  evaluation  of  the  model  performance  we 
compute  the  rms  error  between  the  simulated  SST  and 
the  DMSP  MCSST  (Figure  14).  The  rms  errors  from 
November  13  to  17,  when  cloud  cover  was  diminished, 
ranged  between  0.7°  and  1.3°C.  This  is  little  larger  than 
the  differences  (0.6°C)  between  in  situ  observations  and 
advanced  very  high  resolution  radiometer  (AVHRR)- 
derived  SSTs  [Shay  et  al,  1992], 

6.2.  Maximum  Current  Speed  in  the  Wake 

Previous  studies  [e.g.,  Chang  and  Anthes,  1978;  Great- 
batch,  1983]  show  the  linkage  between  the  storm  trans- 
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Figure  11.  Zonal  cross  section  of  simulated  temperature  anomaly  (degree  Celsius)  relative 
November  4,  1996,  along  12.3°N:  (a)  on  November  16  and  (b)  on  November  17. 


-200 


113  114 


V  (cm/s)  Nov.  9  Local:  13.02  N 


Longitude  (E) 

W  (cm/hr)  Nov.  9  Local:  13.02  N 


Figure  12.  Zonal  cross  section  of  simulated  (a)  v  (cm  s-1)  and  (b)  w  (cm  hr"1)  along  13°N  on 
November  9,  1996. 


lation  speed  and  the  maximum  current  speed  U  in  the 
wake, 

U  =  At  Vt ,  (g) 

where  Vt  is  the  storm  translation  speed  and  AT  mea¬ 
sures  the  ratio  of  the  along-track  advection  terms  to 
the  local  acceleration  terms.  Greatbatch’s  [1983]  esti¬ 
mation  was  higher  than  Chang  and  Anthes’  [1978]  es¬ 
timation  for  the  same  storm  translation  speed.  The 
translation  speed  of  Ernie  is  ~5  m  s'1;  the  maximum 


current  speed  in  the  wake  should  be  2.7  m  s"1  accord¬ 
ing  to  Greatbatch’s  [1983]  estimation  and  '1.7  m  s-1 
according  to  Chang  and  Anthes’  [1978]  estimation.  The 
POM-modeled  maximum  current  speed  was  2  m  s_1, 
which  is  closer  to  Chang  and  Anthes’  [1978]  estimation! 

7.  Conclusions 

1.  This  study  used  the  POM  to  investigate  the  SCS 
responses  to  tropical  cyclone  Ernie  during  November  8- 
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Figure  14.  Root-mean-square  (rms)  errors  (degree 
Celsius)  between  the  simulated  SST  and  the  DMSP 
MCSST. 


18,  1996.  A  tropical  cyclone  wind  profile  model  was 
used  to  simulate  Tropical  Cyclone  Ernie’s  wind  stress 
and  to  force  POM  during  this  14  day  period.  POM  ve¬ 
locity,  temperature,  salinity,  and  surface  elevation  fields 
were  then  analyzed  to  investigate  ocean  responses  pro¬ 
duced  by  the  storm.  Comparison  between  the  modeled 
and  observed  temperature  and  surface  elevation  shows 
the  POM  capability  in  simulating  the  regional  sea  re¬ 
sponse  to  a  moving  tropical  cyclone. 

2.  The  simulation  shows  the  similarities  in  oceanic 
response  to  tropical  cyclones  between  the  western  Pa¬ 
cific  regional  seas  (e.g.,  the  SCS)  and  the  western  At¬ 
lantic/eastern  Pacific  regional  seas:  both  have  strong 
near-inertial,  anticyclonic  turning  upper  ocean  currents 
to  the  right  of  the  storm  track.  These  highly  divergent 
upper  layer  currents  also  generated  the  typical  bias  of 
maximum  SST  cooling  to  the  right  of  the  storm  track. 
The  highly  divergent  surface  currents  produced  strong 
up  welling.  Maximum  cooling  associated  with  this  up- 
welling  (part  of  spreading  three-dimensional  wave  wake) 
was  also  found  to  the  right  of  the  storm  track  at  the  base 
of  the  mixed  layer.  The  model  fields  depicted  intense 
current  shear  between  the  mixed  layer  and  the  thermo- 
cline,  with  180°  reversal  of  these  currents.  Sea  surface 
depressions  developed  in  the  wake  of  the  storm  were 
also  similar  to  studies  of  other  storms. 

3.  The  simulation  shows  several  unique  features  of 
the  SCS  response  to  tropical  storm  Ernie  1996,  such  as 
Kuroshio  intrusion,  coastal  interaction,  and  concentric 
response  during  the  looping  around  stage.  The  POM 
model  successfully  simulated  the  SST  warming  across 
the  Luzon  Strait  from  southwest  Taiwan  to  northwest 
Luzon,  which  was  generated  by  the  convergence  be¬ 
tween  a  northward  flowing  coastal  current  and  the  in¬ 
truded  Kurishio  Current  through  the  Luzon  Strait.  Coas¬ 
tal  interactions  under  storm  forcing  produced  a  subsur¬ 
face  jet  flowing  parallel  to  the  coast  of  Luzon  at  50-100 
m  depth. 

During  Ernie’s  looping  around  stage,  SST  and  sea 
surface  elevation  anomaly  patterns,  which  had  been 
concentrated  in  the  storm’s  wake,  quickly  became  nearly 
concentric  with  the  storm’s  position  and  weakened  as  it 


went  quasi-stationary.  Maximum  SST  cooling,  1®C,  and 
sea  surface  depression,  10  cm,  occurred  near  the  center 
of  the  storm. 

4.  Although  POM  adequately  simulated  the  SCS  re¬ 
sponse  to  Ernie,  much  more  work  in  modeling  ocean 
forcing  by  tropical  cyclones  is  required.  Three-dimen¬ 
sional  observations  of  ocean  temperature,  salinity,  and 
velocity  fields  would  allow  initialization  of  the  model 
with  a  more  realistic  prestorm  thermodynamic  struc¬ 
ture.  Profile  measurements  acquired  during  the  pas¬ 
sage  of  a  tropical  cyclone  would  significantly  improve 
the  analysis  of  the  ocean  responses  to  the  storm  forcing. 
Storm  forcing  would  also  be  more  realistic  by  utilizing  a 
coupled  ocean-atmosphere  model,  such  as  the  Coupled 
Atmosphere  and  Ocean  Coastal  System  (CAOCS)  un¬ 
der  development  by  the  Naval  Postgraduate  School  [Chu 
et  ai,  1999a.]  Use  of  such  a  system  would  provide  more 
accurate  wind  stress  forcing,  through  modifications  of 
the  wind  field  by  surface  frictional  effects,  and  improved 
ocean/atmosphere  thermal  and  salinity  fluxes. 
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